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GEOMETRIC SHAPE FACTORS FOR PLANETARY -THERMAL AND 

PLANETARY-REFLECTED RADIATION INCIDENT UPON 

SPINNING AND NONSPINNING SPACECRAFT 

By Lenwood G. Clark and E. Clay Anderson 
Langley Research Center 

SUMMARY 


Truncated series expressions are obtained for the planetary-thermal and 
planetary-reflected solar radiation incident upon arbitrarily oriented spinning 
and nonspinning spacecraft. For the spinning case, spacecraft having the angu- 
lar momentum vector either normal to or coincident with the axis of symmetry 
are considered. Results are presented in graphical form for planar surfaces 
and bodies of revolution as a function of spacecraft attitude and altitude. 
Results for the spinning case are obtained by averaging instantaneous values 
over a spin period and are sufficiently accurate for most spinning spacecraft. 
Expressions for obtaining the spacecraft’s attitude in terms of orbit param- 
eters for the solar -oriented case are included. 


INTRODUCTION 


In order to maintain the temperatures of many components of spacecraft 
within specified limits, the thermal design of the spacecraft must be based on 
an accurate description of the radiant heat flux incident to the vehicle. The 
sources of radiant heat flux incident upon an orbiting spacecraft are: direct 

solar radiation, planetary-reflected solar radiation, and planetary-thermal 
radiation. The magnitudes of these heat fluxes are influenced by the appropri- 
ate geometric shape factors which vary with altitude and attitude. 

Since an orbiting spacecraft is located at a great distance from the sun, 
the geometric shape factor for direct solar radiation is the easiest to calcu- 
late and is defined by projection of the vehicle area element in question in 
the direction of the sun. However, the geometric shape factors for planetary- 
thermal and planetary-reflected solar radiation are dependent upon both alti- 
tude and attitude and thus are more difficult to determine. 

Geometric shape factors for planetary-thermal and planetary-reflected 
solar radiation are given as integral expressions in references 1 to 4 for 
several geometric configurations. The solution of these equations, however, 
requires the use of digital computers. References 1 to 3 consider various geo- 
metric configurations for both the spin and nonspin modes of operation but 



present a limited amount of data. Reference 4 presents data for the nonspin 
case only. 

This report presents truncated series representations for the geometric 
shape factors for planetary -thermal and planetary-reflected solar radiation 
incident upon arbitrarily oriented planar surfaces and bodies of revolution. 
Each configuration is considered for both the spin and nonspin modes of opera- 
tion and results are presented in graphical form as a function of the space- 
craft's altitude and attitude. The results of this study are of value in the 
thermal design of any orbiting spacecraft, spinning or nonspinning, since they 
can be used to predict the approximate radiant heat flux incident to the vehi- 
cle and thus assist in the design of temperature control systems. 

Expressions for obtaining the spacecraft's attitude in terms of orbit 
parameters for the solar-oriented case are derived in an appendix. 


SYMBOLS 


A surface area 

a planetary albedo, dimensionless 

C s solar constant 

F geometric shape factor, dimensionless 

H distance between spacecraft and center of planet, R + h 

H spacecraft angular momentum vector 

h orbit altitude 

I orbit inclination measured from planet equatorial plane 

i ecliptic inclination measured from planet equatorial plane 

J radiant heat intensity 

M rotational matrix 

N integer 

N unit normal vector to spacecraft surface element 

q radiant heat flux incident upon spacecraft surface area 

R mean radius of planet 

R vector from spacecraft to center of planet 
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a angle "between line from spacecraft to center of planet and line from 

spacecraft to element on planet surface 

(3 right ascension of ascending orbit node measured in planet equatorial 

plane (zero at vernal equinox) 

7 spacecraft position angle measured in orbit plane (zero at ascending 

orbit node), 7 p + x\ 

7 p argument of orbit perigee measured in orbit plane (zero at ascending 

orbit node) 

5 semiapex angle 

€ orbit eccentricity, dimensionless 

spacecraft position angle measured in orbit plane (zero at orbit 
perigee) 

r\ Q angle between line from center of planet to sun and line from center 

of planet to spacecraft 

0 angle between spacecraft axis of symmetry and line from spacecraft 

to center of planet 

A angle between vector normal to spacecraft surface element and line 

from spacecraft to element on planet surface 

A angle between vector normal to spacecraft surface element and line 

from spacecraft to center of planet 

v angle of rotation about either spacecraft orientation vector or 

spacecraft angular momentum vector 

| angle of rotation of spacecraft surface element about spacecraft axis 

of symmetry measured from either spacecraft orientation vector or 
spacecraft angular momentum vector 

p distance between spacecraft and element on planet surface 

0 angle between spacecraft angular momentum vector and line from space- 

craft to center of planet 

cp azimuth angle in spherical coordinates 

£3 position angle of sun measured in the ecliptic plane (zero at vernal 

equinox ) 


Subscripts : 

ent shadow entrance 
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exit shadow exit 

gen general 

L limit 

p perigee 

R planetary-reflected solar 

s solar 

T planetary thermal 


ANALYSIS 


Planetary -Thermal Radiation 

The planetary-thermal radiant heat energy incident upon a vehicle surface 
area A is given by 


q T = JrpA F (l) 

where Jrp is the thermal radiation intensity of a spherical planet assumed to 
radiate according to Lambert’s cosine law. With this assumption 


J 


T “ 


C s (l - a) 

4 


( 2 ) 


Planetary-Re flee ted Solar Radiation 

If it is assumed that the planet is a diffuse reflector of solar radiation 
and that the geometric shape factor for planetary-reflected solar radiation is, 
from reference 5* approximately given by 


F R = F cos t[ b (3) 

then the planetary-reflected ?olar radiant heat energy incident upon a vehicle 
surface area is given by 


q.R = JrAF cos T) s 


W 


k 



where 


Jr — Cga 


J R = ° 



(5) 


The evaluation of equations (l) and (4) requires Imowledge of the geomet- 
ric shape factor for planetary-thermal radiation F. The following section of 
this report is concerned with the derivation of this parameter for several 
spacecraft configurations. 


Geometric Shape Factor 

Planar configuration .- For two black surfaces separated by a nonabsorbing 
medium, the geometric shape factor Fq 2 from surface Aq to surface A 2 is 
defined as the fraction of the total radiant flux leaving Aq that is incident 
upon A 2 and is given for Aq «< A 2 from reference 6 by 


12 = T7 


Ap(A)cos a 


A 2 


cLAp 


7tp* 


( 6 ) 


where Ap(A) is the projected area of surface Aq in the direction of A 2 
and is 


Ap(A) = Axf(A) 


where to a viewer at cLA.2, the area distribution f(A) of surface Ax is 
defined in the interval -it to Jt by 


f (A) = 0 



< A < - 



f(A) = cos A 



A < 



( 7 ) 


f (a) = 0 



< A < 
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By use of the spherical coordinate system of figure 1, equation (6) can be 
expressed as 


f 12 


1 

it 



f(A)sin a da dip 


As in reference 7 , f(A) can be expanded in a Fourier series: 


( 8 ) 


00 



N=0 


cos NA 


(9) 


The series is rapidly convergent and represents f(A) to a high degree of 
accuracy (within 3 percent) when truncated at N = 6. Thus, 


f(A) = — + — cos A + — cos 2A 
it 2 3« 


O p 

— cos 4A + cos 6A 

15n 35* 


( 10 ) 


or equivalently 


f(A) = — + — cos A + — cos^A - cos^A + cos^A 
7* 2 7ir 21n 35 * 


(id 


From figure 1, spherical trigonometry gives the relation 


cos A = cos A cos a + sin A sin a cos cp (12) 

Substitution of equations (ll) and (12) into equation (8) and performing the 
indicated integration gives 

F 12(^> a L) = Bq + Bp cos A + B 2 cos^A + B^ cos^A + B^ cos^A ( 13 ) 
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where 


>v 

P / 57*7 k- 2 )| \ 

330 = 7 ^(l 05 - T cos “L + 3 cos3a L - - COS^ + - cosTojJ 


B-i = ~ sircar 


; 2 = ~^cos a,p - 2 cos^a^ + if cos^ap - 3 cos^a^j 

3 = —• |~cos cos-^ap - ^ cos^a^ + l8 cos^a^j 

4 = 35~ C ° S ” ^5 cos^a L + 63 cos^ap - 33 cos^a^ 


B 


= sln " 1 (| 


H 


= (hp + R)(l + e ) 
1 + € cos rj 


h = H - R 




(14) 


J 


Equation ( 13 ) expresses the instantaneous geometric shape factor between a 
planar surface and a spherical planet as a function of the attitude angle A 
which is defined for an arbitrary orientation in figure 1, and orbit alti- 
tude h. Results of equation ( 13 ) are shown graphically in figure 2. The 
parameters Bq, Bq, Bg, Bj, and B 4 are tabulated for h/R ratios of 0.001 

to 10.0 in table I. For a solar-oriented vehicle, A in terms of orbit param- 
eters (see appendix for derivation) is defined as 


cos A 



( 15 ) 


7 



where 


Ap = cos ft cos p + cos i sin ft sin p 

A2 = cos v sin p sin i + sin v(sin p cos i cos ft - sin ft cos p) 

A^ = cos v(sin p cos i cos ft - sin ft cos p) - sin v sin p sin i 

A ij, = -cos ft cos I sin p + cos i sin ft cos I cos p + sin i sin ft sin I 

A<p = cos v(cos I cos p sin i - sin I cos i) + sin v(sin ft cos I sin p 
+ sin I sin i cos ft + cos I cos p cos i cos ft) 

A^ = cos v(sin ft cos I sin p + cos I cos p cos i cos ft + sin I sin i cos ft) 

+ sin v(sin I cos i - cos I cos p sin i) 


Bodies of revolution .- Geometric shape factors are determined for instan- 
taneous (nonspinning) and spinning bodies of revolution. 


Instantaneous geometric shape factor: If one considers a convex surface 

to be composed of an infinite number of planar elements , the instantaneous 
geometric shape factor for a right circular conical frustum is given by 


F ( 8 ,e,a L ) = i 



where, from figure spherical trigonometry gives 


(16) 


cos A = sin 5 cos 0 + cos 6 sin 0 cos | 


(17) 


Substituting equations ( 13 ) and (17) into equation (l6) gives upon integration 
F( 8 ,e,a L ) = Bq + Bj_(C) + Bg^C 2 + + B 3 (c k + JC 2 !) 2 + | 


+ Bk(c6 + C^D 2 + C 2 D 1+ + -L 

4 \ 2 8 16 / 


(18) 


where 


C = sin 5 cos 0 


D = cos 5 sin 0 
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Equation (l8) expresses the instantaneous geometric shape factor between the 
cone and a spherical planet as a function of the semiapex angle 5, the atti- 
tude angle 0 which is defined for an arbitrary orientation in figure 3, and 
orbit altitude h. The attitude angle 9 in terms of orbit parameters for a 
solar-oriented vehicle is defined in the following equation (derived in the 
appendix): 


cos 0 = 


- jcos /[cos v ( cos i cos ft sin 0 - sin ft cos (3) - sin v(sin i sin pj] 

+ sin 7 [cos v ( sin ft cos I sin (3 + cos i cos ft cos I cos p 

+ sin i cos ft sin i) + sin v(cos i sin I - sin i cos I cos p >3} (19) 


Results of equation (l8) are shown graphically in figures k, 5, and 6. Setting 
5 - 0° or 90° in equation (l8) gives the geometric shape factor for a cylinder 
and a planar surface, respectively. For 5 = 0°, the geometric shape factor 
for a cylindrical configuration is 


= Bo + 


B 


sin20 + 

Q 


sin^0 + iLii sin^0 


( 20 ) 


For 5 = 90°, equation (l8) reduces to equation (13) with 0 = A. Equation (l8) 
also gives the geometric shape for a spinning planar surface where the axis of 
symmetry is the angular momentum vector and the angle 90° - 8 becomes the 
attitude angle of the planar surface with respect to this vector. 


For a general convex body of revolution, the instantaneous geometric shape 
factor is given by 


F(9,a L ) 


L F ( A ’ a B) 


dA 


gen 


( 21 ) 


dA 


where F(x,a L ) is the instantaneous geometric shape factor to an infinitesimal 

planar surface element and dA is the surface element which must be appropri- 
ately described for the general configuration. 

It follows that the instantaneous geometric shape factor for a hemispheri- 
cal configuration can be expressed as 


1 r rt / 2 r 

F(e,a L ) = ~ J I F^A,a£^cos 5 d£ d5 


( 22 ) 
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or equivalently 


ott/2 

F ^0,a L ) = J F^5^9,a L jcos 5 


Substituting equation (l8) into equation (23) and. integrating gives 


F( 0 ,Ol) = Bq+ — cos 9+-— + — + — 


V W ^ 2 3 5 7 

Results of equation (24) are shown graphically in figure 7* 

The geometric shape factor for a spherical configuration is obtained by 
setting 0 = 90° in equation (24) which gives 


(90,CIl) = Bq 


h + h + 


Results of equation (25) are shown graphically in figure 8. 

Spinning geometric shape factor: If the spacecraft is spinning with the 

angular momentum vector normal to the principal axis as shown in figure 9> then 

cos 0 = sin 0 cos v (26) 

If it is assumed that the spin rate is sufficiently fast (that is, the ratio of 
orbit to spacecraft angular velocity is small), the geometric shape factor for 
a cone averaged over a spin cycle is defined as 


f ( 5 ’H = h C F ( 6 ' e ^b 


Substituting equations (l8) and (26) into equation (27) gives upon integration 


F(6,$,a L ) = Bq + ^E 2 + G 2 + 


+ + 4G 2 E 2 + K 2 E 2 + 3G 2 K 2 + 

+ :/e 6 + 9G 2 E^ + 3. eSc 2 + 9E 2 G^ 


+ ZL g 2 k 2 e 2 + 2 e 2 k^ + g^ + gSc 2 + g 2 k^ + 2L 
2 8 2 8 16 
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where 


E = cos 5 
G = sin 6 sin 0 


K 2 = -cos 2 S sin 2 *!* 

Equation (28) expresses the geometric shape factor between a spinning cone and. 
a spherical planet as a function of the semiapex angle 5, the attitude angle $ 
which is defined for an arbitrary orientation in figure 9* and the orbit alti- 
tude h. Results of equation (28) are shown graphically in figures 10, 11, 
and 12. Setting 5=0° gives 



^lLs2$ + 1 
8 l 8 



1 - — sin 2 0 + 2 sin^$ - - 2 - sin^Oj (29) 

v 2 8 16 j 

which is the geometric shape factor for a cylindrical configuration. The atti- 
tude angle $ appearing in equations (28) and (29) is defined (and derived in 
the appendix) in terms of orbit parameters for a solar-oriented spacecraft as 



cos 0 = -[cos /(cos 3 cos 2 + sin 3 cos i sin ft) 

+ sin /(cos I cos 3 cos i sin ft - cos I sin 3 cos 2 
+ sin I sin i sin ft 3 (30) 


Example problem .- As an example problem, consider the manned orbital 
research laboratory concept as shown in the spinning and nonspinning modes of 
operation in figures 13 and l4, respectively. It is required to calculate the 
external radiant heat flux incident upon several points on the surface of this 
vehicle as it orbits the earth for both modes of operation. Also, let the 
vehicle be solar oriented (that is, the axis of symmetry is normal to an orien- 
tation vector which is pointed at the sun) and let the sun and ascending orbit 
node be located at the vernal equinox and let the orbit plane be coincident 
with the ecliptic plane. Thus, for a circular 200-nautical-mile orbit 


0° 

ii 

0 

0 

0° 

C s = 443 Btu/hr-ft 2 

23. 5° 

R = 3^35*058 nautical miles 

23. 5° 

a = O.36 

0 

h = 200 nautical miles 
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The points of interest are the locations of antennae on the cylindrical section 
of the vehicle and the antennae are positioned with respect to the orientation 
vector (pointed at the sun) by £ = 45 ° and £ = 315 ° • (See fig. 3 *) 

Instantaneous (nonspinning) shape factor: The instantaneous geometric 

shape factor is given by equation ( 13 ) where the attitude angle A is defined 
by equation (15). For this mode of operation, assume the axis of symmetry of 
the vehicle to be inclined normal to the ecliptic plane (that is, v = 90 °) > 
and setting 5 = 0 ° (for a cylinder) gives 


cos A = sin | sin 7 - cos | cos 7 


For this case, equation (30) reduces to 


cos 0 = -cos 7 = -cos T) 

‘s 

Substitution of this equation into equations (l) and ( 4 ) gives the planetary- 
thermal and planetary-ref lected solar radiation incident to £ = 45 ° and 
i = 3 b 5 0 for the nonspinning mode of operation shown in figures 15 and 1 6, 
respectively. 

Spinning shape factor: The geometric shape factor is given by equa- 

tion (l8) where 6 = 90° " I and where for this case, the attitude angle 9 
defined by equation (19) is 


9 = l80° - T) g = 180 0 - 7 

Substitution of this equation into equations (l) and ( 4 ) gives the planetary- 
thermal and planetary-ref lected solar radiation incident to | = 43° and 
| = 315 0 for the spinning mode of operation as shown in figures 15 and 16, 
respectively. 

Solar radiation: For | | J < 90 ° , incident solar radiation is given by 


q s = C S A cos | 


in sunlight and by 


q s = 0 


in shadow. Shadow entrance and exit occur when, from reference 7 , 
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cos Tjg = -cos 


where from equation (l4) 


a T = sin-lf— \ = sin~lf — - — 

L \E R + h, 


Solution of these equations gives 

*1 


s ent 


109 • l c 


Vxit = 25 °' 9 

For the example, the incident solar radiation to £ = 45° and £ = 315° 
both modes of operation is 


2s. = 41+3(0.70711) = 313.25 Btu/hr-ft 2 (-109.1° < 7 < 109. 1°) 


<ls 


= 0 


( 109 . i° <; 7 < 250 . 9 0 ) 


where 7 = r) ( 


RESULTS AND DISCUSSION 


The geometric shape factors calculated by the approximate expressions 
given in this report are compared with graphical and tabulated data given in 
references 3, 5, and 8 for several spacecraft configurations in table II. The 
accuracy of the approximate solutions, as demonstrated in table II, is suffi- 
cient for most engineering applications because of the uncertainty regarding 
the variation of local planet surface optical properties from the assumed mean 
values . 

Graphical results for the nonspinning planar surface (eq. ( 13 ) ) are given 
in figure 2. The geometric shape factor for a spinning planar surface is 
determined from equation (l 8 ) and is plotted in figures 4, 5 y an l 6 for non- 
spinning cones with the appropriate semiapex angle 8 . From equation ( 28 ), the 
results for spinning cones having the angular momentum vector normal to the 
axis of symmetry are shown graphically in figures 10, 11, and 12. The geomet- 
ric shape factors for cylindrical bodies are found by setting 8 = 0 ° in equa- 
tions (l 8 ) and ( 28 ). Results for hemispherical configurations (eq. (24)) and 
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spherical configurations (eq. ( 25 )) are shown graphically in figures 7 and. 8, 
respectively. 

Equations (l) and (4) give the planetary -thermal and planetary-ref lected 
solar radiant heat energy incident upon an orbiting spacecraft. If the space- 
craft has solar orientation, the incident radiant heat energy at any point in 
orbit can be found for any desired time after orbit injection with the aid of 
the orbit geometry and the equations presented in the appendix. 

The radiant heat energy incident upon configurations not presented in this 
report may be found by approximating the vehicle surface by planar elements. 

If it is desired to determine the peripheral variation of radiant heat 
energy upon the configurations presented in this report, the appropriate geo- 
metric shape factor expression can be integrated between the desired peripheral 
positions, or the configuration may be approximated by planar surface elements. 


CONCLUDING REMARKS 


Truncated series expressions for the geometric shape factors for planetary- 
thermal and planetary-reflected solar radiation incident upon arbitrarily 
oriented spinning and nonspinning spacecraft are derived and discussed. 

Results indicate good agreement between the approximate closed form solutions 
of this report and existing computer solutions found in the literature. 

The results of this study are of value in the thermal design of any 
orbiting spacecraft, spinning or nonspinning, since they can be used to predict 
the approximate radiant heat flux incident to the vehicle without the use of 
computers, and thus assist in the design of temperature control systems. 

The radiant heat energy incident upon configurations not presented in this 
report may be found by approximating the vehicle surface by planar elements. 


Langley Research Center, 

National Aeronautics and Space 
Langley Station, Hampton, 


Admini strati on , 

Va., March 2, 1965 . 
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n terms of orbit parameters 
>y considering a nonrotating 
,1 planet with a fixed right 
.hogonal coordinate system 
id Z, with the Z-axis alined 
t planet’s axis of rotation 
’ole) and the X-axis alined 
■ intersections of the equa- 
nd ecliptic planes (vernal 
. (See sketch (a).) 

1 position of a spacecraft 
to the fixed coordinate 
s determined by considering 
nate system X’, Y’, and Z‘ : 
ie of the orbit having the 
perpendicular to the orbit 
id the X’-axis alined with t. 
.ft’s position as shown in 



APPENDIX 


Let i, j , and k represent unit vectors in the X, Y, and Z directions 
and i ’ , y, and k‘ unit vectors in the X', Y' , and Z' directions. The unit 
vectors i 1 , and k' are defined by 


where 


R 


R 

y 

> = M (7) M(I) M(P)< 

J 

k 1 

v. 


k 

V- J 


M(/) 


cos 7 sin 7 0 

-sin 7 cos 7 0 

0 0 1 


M(I) 


M(p) 


10 0 


0 cos 

I 

sin 

I 

0 -sin 

I 

cos 

I 

- 



- 

cos p 

sin 

P 

— 

0 

-sin p 

cos 

P 

0 


0 0 1 


Matrix multiplication gives 



cos 7 cos f 3 - sin 7 cos I sin p 
-sin 7 cos g - cos 7 cos I sin p 
sin I sin (3 


cos 7 sin 3 + sin 7 cos I cos 3 

sin 7 sin I 



sin 7 sin p + cos 7 cos I cos p 

cos 7 sin I 

< 

1 

-sin I cos p 

cos I 


L k J 


The position of the sun relative to the fixed coordinate system is deter- 
mined "by considering a coordinate system X " , Y" , and Z" in the plane of the 
ecliptic having the Z"-axis perpendicular to the ecliptic plane and the X"-axis 
alined with the sun as shown in sketch (c). 
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APPENDIX 



Unit vectors in the X " , Y " , and Z" coordinate system are given by 


where 



V 



n 

i 


< 

j” I 

= M(ft) M(i)< 

j 

> 


k" 

L J 



k 




cos ft 

sin 

ft 

0 

M(n) = 

-sin ft 

cos 

a 

0 



0 

0 


1 



1 0 


0 

-i 

M(i) = 

0 cos 

i 

sin 

i 



0 -sin 

i 

cos 

i 

_ 


IT 


APPENDIX 


Therefore , 


r 


r~ 









r ^ 

1 1 


cos 

n 

cos i 

sin 

Q 

sin i 

sin 

a 


1 

j" 

> ■ 

-sin 


cos i 

cos 

Q 

sin i 

cos 

a 


j 

k” 


0 


-sin 

i 


cos 

i 



k 

L J 


_ 









L J 


If the spacecraft is rotated about the solar unit vector i" through 
the angle v, the unit vectors in a sun-oriented vehicle coordinate system 
X’ ' 1 , Y* 1 ’ , and Z T ‘ ' are defined by 








L k "j 


l k "J 


where 


10 0 
M(v ) = 0 cos v sin v 

0 -sin v cos v 


Then i 11 ’, j 1 1 ’ and k* * ’ are defined in terms of the fixed coordinate system 
by the transformation 


r ^ 


r 

i T ' ' 


i 

J’" 

) = M(v) M(fi) M(i)<( 

j 

k' 1 ' 


k 

v J 


L J 


whereby matrix multiplication gives 


fi 



... 


cos ft 



cos 

i 

sin ft 


sin 

i sin ft 




r 

i 

-cos v sin 

ft 

COS V 

cos i 

cos 

ft - sin i 

sin v 

cos v sin i 

cos ft + sin v 

cos 

i 

< 

j 

sin v sin 

ft 

-sin v 

cos i 

cos 

ft - cos V 

sin i 

-sin v sin i 

cos ft + cos V 

cos 

i 


It 


18 



APPENDIX 


Let e^, e^} and be unit vectors in the X^_, Xg, and X^ directions 

of the right-hand orthogonal vehicle coordinate system shown in figure 3, and 
then the unit vector normal to a planar surface element in this system is given 

by 


N = ep cos 5 cos £ + e 2 sin 8 + e-z, cos 8 sin | 


For a solar-oriented vehicle, if e^_ is selected as the orientation unit vec- 
tor and e 2 as the axis of symmetry for the vehicle configuration, then 


ep = i 


e 2 = J * ' * 


e ? = k' * ' 


and the normal vector becomes 


N = i ,,! cos 8 cos | + j ,,r sin 8 - k f 1 1 cos 8 sin | 


(Al) 


For an arbitrary solar-oriented spacecraft, the attitude angles A, 0, 
and $ can then be defined by 


cos 

>* 

II 

-(N • 

i’) 

cos 

e = ■ 

-(j 11 

' • i ' ) 

cos 

$ = 

“( i 1 T 

' . i«) 


Performing the indicated mathematical operation yields 
cos A 


= - jj 2 os cos | - A 2 sin |)cos 8 + A^ sin sj 

+ sin A^A^ cos | - Atj sin |) cos 8 + Ag sin &jj 


(A2) 
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APPENDIX 


where 


A]_ = cos ft cos p + cos i sin ft sin p 

A2 = cos v sin p sin i + sin v(sin p cos i cos ft - sin ft cos p) 

A^ = cos v(sin p cos i cos ft - sin ft cos p) - sin v sin p sin i 

A^ = -cos ft cos I sin p + cos i sin ft cos I cos p + sin i sin ft sin I 

A^ = cos v(cos I cos p sin i - sin I cos i) + sin v(sin ft cos I sin p 
+ sin I sin i cos ft + cos I cos p cos i cos ft) 

A^ = cos v(sin ft cos I sin p + cos I cos p cos i cos ft + sin I sin i cos ft) 

+ sin v(sin I cos i - cos I cos p sin i) 

and 


cos 0 = 


- ^cos 7 [cos v(cos i cos ft sin p - sin ft cos p) - sin v(sin i sin p)] 
+ sin 7 [cos v(sin ft cos I sin p + cos i cos ft cos I cos p 
+ sin i cos ft sin I) + sin v(cos i sin I - sin i cos I cos 3 jj} 


cos 0 = - [cos 7(cos p cos ft + sin p cos i sin ft) 

+ sin 7(cos I cos p cos i sin ft - cos I sin p cos ft 
+ sin I sin i sin ft 0 (AM 

For an elliptical orbit, the spacecraft position angle 7 (measured from 
ascending orbit node) can be replaced for convenience in the previous attitude 
angle equations by the angle 7 p + where 7 p is the argument of the perigee 

(measured from ascending orbit node) and q is the spacecraft position angle 
measured from the orbit perigee. The orbit altitude, for this case, is given 
by equation (lb) as a function of the perigee altitude h p and orbit eccen- 
tricity e. With this substitution cos $ becomes - cos r) s . 
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TABLE I.- PARAMETERS DEFINED BY EQUATION (l4) TABULATED 
FOR h/R RATIOS FROM O.OOl TO 10.0 


h/R 

Bo 

Bi 

b 2 

b 3 

*4 

0.001 

0.47133 

0 . 49900 

0 . 01619 

-0.00791 

0.01603 

.002 

. 45960 

. 49801 

.02279 

-.01088 

.02234 

.003 

. 45062 

. 49701 

.02778 

-.01295 

.02695 

.004 

. 44308 

. 49602 

.03193 

-.01453 

.03066 

.005 

•43646 

.49504 

.03554 

-.01578 

.03377 

.006 

. 43050 

. 49405 

.03875 

-.01678 

.03644 

.007 

. 42503 

.49307 

.o4i66 

-.01758 

.03877 

.008 

. 41996 

. 49210 

.04434 

-.01822 

. 04082 

.009 

. 41521 

.49112 

.04682 

-.01873 

.04264 

.01 

.41074 

.49015 

.04913 

-.01912 

.04427 

.02 

.37524 

.48058 

.06654 

-.01889 

•05353 

.03 

.34888 

. 47130 

.07825 

-.01443 

.05565 

.04 

.32737 

.46228 

.08699 

-.00794 

.05403 

.05 

. 30901 

.45351 

. 09386 

-.00045 

.05023 

.06 

.29292 

.44500 

.09947 

.00746 

.04508 

.07 

.27858 

.43672 

. 104l6 

.01541 

.03915 

.08 

.26562 

. 42867 

. 10819 

.02316 

.03279 

.09 

.25382 

. 42084 

. 11170 

.03054 

.02626 

.10 

.24298 

.41322 

.11483 

.03747 

.01972 

.20 

. 16729 

.34722 

. 13609 

.07596 

-.03032 

.30 

.12301 

.29586 

. 15017 

.07118 

-.04807 

.40 

.09396 

.25510 

.15976 

.04693 

-.04617 

.50 

.07374 

.22222 

.16514 

.01785 

-.03571 

.60 

.05912 

.19531 

. 16691 

-.00934 

-.02279 

.70 

.04824 

.17301 

.16583 

-.03220 

-.01029 

.80 

.03996 

.15432 

.16268 

-.05026 

.00061 

.90 

.03353 

. I385O 

. 15810 

-.06392 

.00960 

1.0 

.02847 

. 12500 

. 15260 

-.07384 

.01674 

2.0 

.00849 

.05556 

.09456 

-.08116 

•03359 

3.0 

.00394 

.03125 

.05942 

-.05787 

.02588 

4.0 

.00227 

.02000 

.03999 

-.04098 

.01886 

5.0 

.00148 

.OI389 

.02853 

-.03001 

.01400 

6.0 

.00105 

.01020 

.02130 

-.02275 

.01070 

7.0 

.00078 

.00781 

.01648 

-.01777 

.00840 

8.0 

.00060 

.00617 

.01312 

-.01423 

.00675 

9.0 

.00048 

.00500 

.01068 

-.01164 

.00553 

10.0 

.00040 

.00413 

.00886 

-.00969 

.00461 
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TABLE II.- COMPARISON OF RESULTS 


Configuration 

h/R 

Attitude, deg 

F 

(*) 

F 

(**) 

Reference 

Nonspinning 

0.1164 

A = 

0 

0.8043 

0.8020 

8 

planar surface 

.1164 

A = 

60 

. 46 o 4 

.4611 

8 

(eq. (13)) 

.1164 

A = 

90 

•0593 

.0598 

8 


O.2911 

A = 

0 

0.6004 

0.6000 

8 


.2911 

A = 

60 

.3172 

.317*+ 

8 


.2911 

A = 

90 

• 0173 

.0172 

8 

Nonspinning body 

0.1164 

5 = 0 

6=0 

0.2269 

0 . 2267 

5 and 8 

of revolution 

.1164 

5 = 0 

0 = 30 

.2429 

.2419 

3 

or spinning 

.1164 

5 = 0 

0 = 90 

.3071 

.3054 

5 and 8 

planar surface 

.1164 

6 = 0 

0 = 90 

.3071 

.3046 

3 

(eq. (18)) 

.1164 

6 = 60 

0 = 30 

.6219 

.6182 

3 


0 . 2911 

6=0 

0=0 

0 . 1262 

0 . 1262 

5 and 8 


.2911 

6 = 0 

0 = 30 

.1463 

.1456 

3 


.2911 

5 = 0 

0 = 90 

.2131 

.2126 

5 and 8 


.2911 

5 = 0 

0 = 90 

.2131 

.2105 

3 


.2911 

6 = 60 

0 = 30 

.4544 

.4569 

3 

Hemisphere 

0.1164 

0 = 

0 

0 . 4781 

0 . 4768 

5 and 8 

(eq. ( 24 )) 

.1164 

0 = 

60 

.3778 

.3850 

5 and 8 


.1164 

0 = 

120 

.1772 

.i860 

5 and 8 


.1164 

0 = 

l80 

.0769 

.0788 

5 and 8 


0 . 2911 

0 = 

0 

0.3336 

0.333 1 * 

5 and 8 


.2911 

0 = 

60 

.2386 

.2629 

5 and 8 


.2911 

e = 

120 

.1086 

.1134 

5 and 8 


.2911 

0 = 

180 

.0336 

.0344 

5 and 8 

Sphere 

0.1164 



0.2775 

0.2778 

8 

(eq. ( 25 )) 

.2911 



.1836 

.1840 

8 


.5000 



.1272 

.1274 

5 


.7000 

- 



.0956 

.0956 

5 


Values calculated by methods of this report. 
**Values taken from indicated reference. 
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Figure 2.- Instantaneous geometric shape factor for a planar surface as a function of altitude and attitude 
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Figure 4.- Instantaneous geometric shape factor for a right circular cylinder (5 = 0°) as a function of altitude 

and attitude . 




(a) 6 - 10°. 

Figure 5*~ Instantaneous geometric shape factor for a cone as a function of altitude and attitude 
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(b) 8 = 20°. 


Figure 5.- Continued. 



(c) 5 - 30°. 


Figure 5-' Continued 
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(d) 5 = 1+0°. 
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Figure 5 *~ Continued 
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Figure 5 


Continued 
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Figure 5.- Continued 
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Figure 6.- Instantaneous geometric shape factor for a flat plate (5 = 90°) as a function of altitude and attitude 



Figure 7* - Instantaneous geometric shape factor for a hemisphere as a function of altitude and attitude. 
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Figure 8.- Geometric shape factor for a sphere as a function of altitude 




Figure 10.- Geometric shape factor for a spinning right circular cylinder (5 = 0°) as a function of altitude and 

attitude . 
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(a) 5 = 10°. 

Figure 11.- Geometric shape factor for a spinning cone as a function of altitude and attitude. 
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(b) 6 = 20°. 


Figure 11.- Continued. 
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(c) 5 = 30°. 
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Figure 11. ~ Continued 
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(d) 6 = l0°. 

Figure 11.- Continued. 
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(e) 6 = 50 °. 
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(f) s - 6o°. 


Figure 11.- Continued 
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(g) 5 = 70°. 
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Figure 11.- Continued 
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Figure 12.- Geometric shape factor for a spinning flat plate (5 = 90°) as a function of altitude and attitude. 



Figure 13 . - Spinning configuration of a manned orbital research laboratory. 
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Figure 14.- Nonspinning configuration of a manned orbital research laboratory. 
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Figure 15.- Incident earth-thermal radiation to antennae locations as a function of orhit position for spinning and 

nonspinning modes of vehicle operation. 





